Polarized Electric Current in Semiclassical Transport with Spin-Orbit Interaction 
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Semiclassical solutions of two-dimensional Schrodinger equation with spin-orbit interaction and 
smooth potential are considered. In the leading order, spin polarization is in-plane and follows the 
evolution of the electron momentum for a given subband. Out-of-plane spin polarization appears 
as a quantum correction, for which an explicit expression is obtained. We demonstrate how spin- 
polarized currents can be achieved with the help of a barrier or quantum point contact open for 
transmission only in the lower subband. 
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I. INTRODUCTION 

Achieving spin manipulation in nanodevices by means 
of electric fields (without using less selective magnetic 
fields) represents the ultimate goal of spintronics. Spin- 
orbit interaction, which couples electron momentum to 
its spin, is one of the most promising tools for realiz- 
ing spin-polarized transport [lj, [2j . Several schemes lead- 
ing either to spin accumulation or to polarization of the 
transmitted current induced by the spin-orbit interac- 
tion have been put forward. Predictions of electric field 
induced spin-accumulation at the boundaries of a sam- 
ple, which originates from asymmetric scattering from 
impurities [3,0] (extrinsic spin-Hall effect) or from spin- 
orbit split band structure (intrinsic effect), has re- 
cently reached a stage of experimental realization |7j . In- 
plane bulk spin polarization appears in the electric field 
in two-dimensional systems with broken inversion sym- 
metry Spin polarization in quantum wires with low 
carrier density has been shown to occur due to the in- 
terfaces of spin-degenerate and spin-split regions [j| . In- 
terfaces between two-dimensional regions with different 
spin-orbit splitting have also been used for that purpose, 
in the case of a sharp [H EU or an arbitrary Il2| inter- 
face, as was the scattering from a sample edge |l3l IT^. 
Other proposals include polarization due to tunneling 
through a double-barrier structure |l5l IT^ | and tunnel- 
ing between two quantum wires |17| . Reference |l8j sug- 
gested a three-terminal device with a spin-orbit split cen- 
tral region as a spin filter, which was numerically tested 
by Refs. 0] and [2(j ■ Reference |2l| pointed to a pos- 
sibility of generating spin-polarized currents by utilizing 
crossings of spin-orbit-split subbands belonging to differ- 
ent transverse channels. These proposals are still lacking 
experimental realization. 

In the present paper we suggest a way to polarize elec- 
tric currents by passing them through a region where, by 
increasing the external electrostatic potential, the upper 
spin-orbit-split subband is locally positioned above the 
Fermi level. The proposed method utilizes electric gat- 
ing whose effect is two fold: (i) it completely suppresses 
transmission via the upper spin-orbit-split subband, and 



(ii) it allows transmission only in a narrow interval of in- 
cident angles in the lower subband. In contrast to the 
proposals which advocate strong variations of the spin- 
orbit coupling and, thus, rely on strong gate voltages, 
our method requires only weak potentials of the order 
of a few millivolts (which is a typical scale of the Fermi 
energy). In addition, we predict a specific pinch-off be- 
havior of the conductance, which would allow to detect 
polarized currents without actual measurement of spin. 

We consider ballistic electron transport in gated two- 
dimensional electron gas with the Hamiltonian 

mA 2 

H = — + X(p y (T x - p x a y ) + +V(x,y). (1) 

For the sake of simplicity we concentrate on the case of 
the "Rashba" spin-orbit interaction (the same method, 
however, can be used for more complicated interactions). 
Construction of semiclassical solutions of the Schrodinger 
equation with the Hamiltonian Q follows the reasoning 
of the conventional WKB approach HHHHHHl, which 
is valid for a smooth potential, fr\ VV| min(p 3 /m,p 2 A). 
The advantages of semiclassics are twofold. First, it al- 
lows us to obtain approximate analytical solutions for 
otherwise complicated problems. Second, as we will see, 
it turns out to be especially simple to achieve strong 
polarization of electron transmission in the semiclassical 
regime. 

The Mexican hat shape of the effective kinetic energy 
in the case of spin-orbit interaction leads to a variety of 
unusual classical trajectories (see Fig. 2 below), which 
have never been investigated before. Our approach em- 
ploys strong spin-orbit interaction (or smooth external 
potential) sufficient to affect individual electron trajec- 
tories, in contrast to previous semiclassical treatments 
1 2(1 127| which consider spin-orbit interaction as a pertur- 
bation. Still we do not require the spin-orbit interaction 
to be comparable with the bulk value of the Fermi energy. 
To produce spin-polarized current, it will be sufficient to 
make spin-orbit interaction comparable with the kinetic 
energy at some particular area of the system, for exam- 
ple, near the pinch-off of a quantum point contact. 
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II. SEMICLASSICAL WAVE FUNCTION 

Without the external potential V, the electron spec- 
trum consists of the two subbands, E±(p x ,p y ) = (p ± 
m\) 2 /2m. The subbands meet at only one point, p = 0, 
and the spin in each subband is always aligned with one 
of the in-plane directions perpendicular to the momen- 
tum p. The semiclassical electron dynamics |22j naturally 
captures the essential features of this translationally in- 
variant limit. The classical motion in each subband is 
determined by the equations of motion which follow from 
the effective Hamiltonian: 

H«=<**£? + V(x,v). (2) 
Ira 

Despite the fact that spin does not appear in this equa- 
tion, one can easily construct semiclassical wave func- 
tions, which have spin pointed within the xy plane per- 
pendicular to the momentum: 

V 2p V ±\/Py ~ IPx J 

Here the action S is related to the momentum by p — 
VS, and p — v)u is the classical density for a family 
of classical trajectories corresponding to a given energy 
E. The action S obeys the classical Hamilton- Jacobi 
equation, (|V5| ± mX) 2 /2m + V — E. Application of 
the Hamiltonian to the approximate wave function 
ipo gives, after some algebra, 

ih 

Hi/jo = £tyo - — (V • pv) + h\Fa z ij , (4) 
2p 

where (summation over repeating indices is assumed) 

The second term in the rhs of Eq. (0J vanishes due to the 
continuity equation 

V • pv — 0, v=p/m±\p/p. (6) 

The last (~ a z ) term in indicates that the spin of 
an accelerated electron cannot exactly stay in the plane 
of propagation and acquires a small ~ hW projection 
onto the z axis. To take into account this out-of-plain 
spin precession one has to go beyond the approximation 
of Eq. ©, which is done by 

V; = (1 + hfa z )ip - (7) 

Since (H — E)fa z ipQ — ^f2\pfa z ipo, to the lowest order 
in h, one can relate the functions F and / 

/ = ±F/2p, (8) 

and find the out-of-plane spin density [F is found from 
Eq. ©] 

= ±^F. (9) 
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FIG. 1: Conductance (in units of e 2 /h), and spin polarization 
of the current vs gate voltage (in units of mA 2 /2). Dashed 
lines show the smoothed curves (1141151 . solid lines show the 
quantized values for mXL/h — 10.57T. Dotted line shows the 
conductance without spin-orbit interaction. 

Note that Eq. does not describe the nonadiabatic 
transitions between subbands. After the electron leaves 
the region with nonzero potential gradient, W ^ 0, the 
in-plain spin orientation is restored. 

The out-of-plane polarization of the electron flow in 
the external potential is a subject of the rapidly develop- 
ing field of the spin-Hall effect [1 II II S, . Our result, 
Eqs. JSJ and 10, incorporates previous calculations of 
Ref. |(| which were restricted to the one-dimensional form 
of the potential, V(x), with p y being the integral of mo- 
tion. The validity of Eq. however, is not restricted to 
a simple one-dimensional case and describes the out-of- 
plain polarization for any smooth two-dimensional poten- 
tial (including confining potentials which create quantum 
wires, quantum dots, etc.). In particular, Eq. @ may 
serve as a good starting point for an analytical calcula- 
tion of the edge spin accumulation in ballistic quantum 
wires j28[ [2^. We leave further investigation of these 
interesting effects for subsequent research. 

Solutions of the form, Eq. J3J|, have clear and im- 
portant consequences. During its motion, an electron 
changes the momentum p but always remains in the same 
spin-subband. To change the subband the electron tra- 
jectory should pass through the degeneracy point where 
both components of momentum vanish simultaneously, 
p = 0, which is generically impossible. Moreover, with 
the proper use of potential barriers, one may realize a 
situation where electrons of only one subband are trans- 
mitted and the others are totally reflected. This leads to 
strong polarization of the transmitted electron flow. 

III. SHARVIN CONDUCTANCE 

To give an example of such a spin-polarized current let 
us consider transmission through a barrier, V{x), varying 
along the direction of a current propagation. We assume 
periodic boundary conditions in the perpendicular direc- 
tion (y + L = y). As such a condition makes p y the 
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integral of motion, mixing of orbital channels, which is 
strongly suppressed for generic smooth potential J3J), is 
now absent exactly. For a smooth potential V(x) the con- 
duction channels may either be perfectly transmitting or 
completely closed. The conserved transverse momentum 
takes the quantized values, p™ = 27rhn/L. Consider the 
functions 

For n / the function E±(p x ) splits into two distinct 
branches. At any point x the equation 

El{p x ) =E F -V(x) (11) 

yields solutions p x and p x , corresponding to left- and 
right-moving electrons. Application of a small bias im- 
plies, e.g., the excess of right movers over left movers far 
to the left from the barrier. Particles are transmitted 
freely above the barrier if Eq. fl"l"|) has a solution, p x , for 
any x. Let /i = Ep — V max be the difference between the 
Fermi energy and the maximum of the potential. The 
nth channel in the upper branch opens when 

H = (2tt%| + mXL) 2 /2mL 2 . (12) 

For the lower branch E™(p x ) Eq. i|ll|) has four solutions 
(two for right and two for left movers) for \n\ < m\L/2Trh 
and x close to the top of the barrier. However, far from 
the barrier (where the excess of right-movers is created) 
there are still only two crossings described by Eq. ifTTjl. 
one for right and one for left movers. As a result, all the 
extra electrons injected at x = — oo follow the evolution 
of a solution of Eq. (|llf) with the largest positive p x . For 
all \n\ < mXL/2nh such a solution does exist for any pos- 
itive /i. Thus, at /i = as many as uq = mXL/Trh chan- 
nels open up simultaneously. The channels with higher 
values \n\ > mXL/2nT% in the lower subband E™ open 
when 

fi= (2-Kh\n\ -mXL) 2 /2mL 2 . (13) 

According to the Landauer formula, ballistic conductance 
is given by the total number of open channels multiplied 
by the conductance quantum Go = e 2 /h 

n _ n L / y/2jlm + m\, 0<fi<mX 2 /2 , 
G ~ G °^\ 2^2^, ^>mX 2 /2. ^ 

This dependence G(/i) is shown in Fig. 1. The strik- 
ing evidence of the presence of spin-orbit interaction is 
the huge jump of the conductance at the pinch-off point, 
as opposed to the conventional square-root increase in 
the absence of spin-orbit coupling. This jump is a con- 
sequence of the "Mexican-hat" shape of the spectrum 
E^(p x ,p y ). Accuracy of Eqs. (|T2|) and IT5|) is sufficient 
to resolve the steps in the conductance due to the discrete 
values of |n| = 0, 1,2, (conductance quantization), as 
shown in Fig. 1. The steps in G(fi) are abrupt in the 
limit dV/dx -> 0. 



Close to the pinch-off, at /i < mX 2 , the conserved p y 
component of the electronic momentum varies for differ- 
ent transmitted channels within the range \p y \ < m\. 
Therefore, far from the barrier, where the Fermi momen- 
tum is large pp 3> mX, we have p x 3> p y and transmit- 
ted electrons propagate in a very narrow angle interval 
\0\ < -J mX 2 /2Ep <C 1. Since the electron spin is perpen- 
dicular to its momentum, we conclude that the current 
due to electrons from each of the subbands is almost fully 
polarized. The total polarization of the transmitted cur- 
rent is given by the difference of two currents 



(a y ) = (^a y v x ip)/(ip^v x ij)) = min(l, y / mX 2 /2fi), (15) 

which is also depicted in Fig. 1. This current polarization 
may also be viewed as a creation of in-plain nonequilib- 
rium spin density, maximal on the barrier. 

Vanishing transmission for electrons from the upper 
band for < /i < mX 2 /2 Ijl4(l resembles the total inter- 
nal reflection suggested for creation of polarized electron 
beams in Ref. [IJ- Unlike the latter case, in our proposal 
there is no need to collimate incident electron flow, since 
the upper band electrons are reflected at any angle. 

Semiclassical formulas (|14|l and (|15fl are valid provided 
that there are many open transmission channels, and 
account correctly for the electrons with ^ 0. The 
case n = 0, however, requires special attention. The 
curve E±(p x ) does not split into the lower and upper 
branches, but instead consists of two crossing parabolas 
shifted horizontally. Right movers from both parabolas 
are transmitted or reflected simultaneously. The electron 
flow due to the channels with n = is, therefore, unpo- 
larized. For small n / the crossing of two parabolas 
is avoided. However, the electrons from the upper sub- 
band E" may tunnel into the lower branch E™ in the 
vicinity of the point p x = 0, which results in the decrease 
of spin-polarization of the current. Let the barrier near 
the top has a form V(x) = —mfl 2 x 2 /2. Simple estima- 
tion shows that classically forbidden transition between 
the subbands do not change the net polarization of the 
current as long as hH, <C mX 2 . 

Our results Eqs. i|14fl and (|15fl were obtained for the 
periodic boundary conditions. However, the boundary 
conditions do not play important role for the conductance 
(G oc L) if the width of the "wire" is large compared with 
the width of the barrier, i.e., if L ^> y/h/mfl h/mX. 
If the transverse confinement in the wide wire is ensured 
by the smooth potential [30j the semiclassical transmit- 
ted scattering states may be constructed explicitly using 
the method of Ref. [34j . However, since the spin-orbit 
interaction in our approach appears already in the clas- 
sical Hamiltonian @, calculation of smoothed conduc- 
tance l|14|l requires only a simple counting of classical 
trajectories |3jJ. Our next example below demonstrates 
such semiclassical treatment of realistic boundary condi- 
tions. 
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IV. QUANTUM POINT CONTACT 

Let us consider probably the most experimentally rel- 
evant example of a quantum point contact, described by 
the potential 



V(x,y) = - 



m£l 2 x 2 mu! 2 y 2 



(16) 



We will see that even in this simple model the electron 
flow in the presence of spin-orbit interaction acquires a 
number of interesting and peculiar features. Classical 
equations of motion follow in the usual manner from the 
effective Hamiltonian (J2J: f = dH e g/dp, p = —dH c g/dr. 
We consider quantum point contact (QPC) close to the 
opening with only the lower E— subband contributing to 
the conductance. A crucial property of the Hamiltonian 
H e g, Eq. (J2J, is the existence of a circle of minima of the 
kinetic energy at \p\ = mX. Expanding around a point on 
this circle, p Xo = mAcosa, p yo — mAsina, one readily 
finds the equations of motion for V = p x cos a+p y sin a — 
mX -C toA, 

V + (-ft 2 cos a 2 +uj 2 sma 2 )V = , a = 0. (17) 

The trajectory is found from the relations, x = 
■Pcosa/ra , y = V sina/m. We observe from Eq. 1171) 
that only the trajectories within the angle 



tan|a| < tanag = Cl/u> 



(18) 



are transmitted through QPC. Trajectories with larger 
angles arc trapped (oscillate) within the point contact. 
Examples of both types of trajectories are presented in 
Fig. 2. Quantization of trapped trajectories would give 
rise to a set of (extremely) narrow resonances in the con- 
ductance, specific for spin-orbit interaction. We leave the 
detailed investigation of these narrow features for future 
research. Below we consider only the smoothed conduc- 
tance. 

To calculate the current J through QPC one has to 
integrate over the phase space of the states which are 
transmitted from left to right, 



J 



dy 



d 2 p 



= GV, 



(19) 



and have the energy within the interval fi — eV/2 < E_ < 
H + eV/2, with V standing for the applied voltage. In this 
section we define /i as the difference between the Fermi 
energy and the value of the potential at the saddle point 
fi = Ef — V(0, 0). The integral is most simply evaluated 
at x = (with the velocity given by v x = Vcosa/m). 
The allowed absolute values of the momentum are 



2(i - eV - muj z y z < V /m < 2[i + eV - rm) y*. (20) 

The angle interval of transmitting trajectories consists of 
two domains: \a\ < olq, V > 0, and \a — w\ < a>o, V < 0. 
The appearance of the latter range of integration is highly 



non-trivial. A simple reasoning shows that the particles 
with the velocity antiparallcl to the momentum (v x > 0, 
p x < 0) should not contribute to the conduction in the 
case of a transition through a one-dimensional barrier 
V = V(x), see Eq. I|14ll . Despite corresponding to the 
right-moving electrons, these states do not originate in 
the left lead. Indeed, they exist only in the vicinity of 
x = 0, but disappear as x — > — oo and, thus, cannot 
be populated by the excess electrons (except due to the 
tunneling transitions which are irrelevant in the semi- 
classical regime). Such trajectories, however, do exist in 
QPC, Eq. ifrSJl. as demonstrated in Fig. 2. After passing 
through QPC the trajectory bounces at the wall reversing 
its velocity. This kind of classical turning points, where 
both components of the velocity vanish simultaneously, 
are specific for the effective Hamiltonian The exis- 
tence of transmitting trajectories with | a— ir\ < oto, g < 
results in the doubling of the conductance. Simple cal- 
culation yields 



G — Go 



4mAsina;o /2/i 
ttTiuj 



(21) 



The presence of a threshold angle ao, as well as the 
square-root dependence of G(fi), are in a sharp contrast 
to the well-known result G = Gafj,/Tr?iuj, in the absence 
of spin-orbit interaction. 

Equation (|21|) is valid in the case of many open chan- 
nels. Since Eq. I|17fl describes only the linearized elec- 
tron dynamics, Eq. (|21|) is formally valid if p, mX 2 . 
Nevertheless, the current remains totally polarized for 
< fi < mX 2 /2 [similar to Eq. JTSJ] 



0~yV X 



(22) 



With increasing the chemical potential, ^ > mX 2 /2, 
transmission via the upper subband E + kicks in and the 
degree of polarization gradually decreases, similarly to 
Eq. (|15f) . though with different, more complicated, de- 
pendence of spin-polarization on fi. Note that transmis- 
sion of different orbital channels through QPC is indepen- 
dent as long as the confining potential (|16fl is smooth over 
a distance of the characteristic spin-orbit length Ti/mX. 
It is easy to see that this requirement is equivalent to the 
condition that (u>, O) <C mX 2 /%. This is also a condition 
of large conductance G>Gq. 



V. DISCUSSION 

In both analyzed systems (of ballistic Sharvin conduc- 
tance and of QPC) polarization of current is achieved 
when many channels arc transmitting. As a consequence 
of the Kramers degeneracy, transmission eigenvalues al- 
ways appear in pairs in the presence of time-reversal sym- 
metry, leading to the prohibition of the spin-current in 
the lowest (n = 0) conducting channels (cf. Ref. [lg|). 
In the case of higher channels, however, the degenerate 
transmission eigenvalues belong to the same spin-orbit 
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FIG. 2: Three kinds of trajectories in the point contact, a, 
transmitted trajectory whose momentum is always collinear 
with the velocity. b, trajectory bouncing inside the QPC. 
This trajectory is periodic in the linearized approximation 
described in the text, while the exact calculation for finite am- 
plitude shows its slow drift, c, transmitted trajectory whose 
momentum inside the contact is opposite to the velocity. Elec- 
trons flow from left to right. Arrows show momentum and 
spin orientations. Few equipotential lines are also shown. 

subband and carry, respectively, the same spin polariza- 
tion. For example, in the case of the QPC any trans- 
mitted trajectory x(t), y(t) (e.g., one of the two shown in 
Fig. 2) is accompanied by its mirror reflection x(t), —y(t) 
with identical transmission. 

In InAs-based heterostructures, typical value of spin- 
orbit coupling [sBl is A7i = 2 x 10~ 11 eVm. Character- 



istic spin-orbit length Ir = Ti/m*\ = 100 nm and en- 
ergy m*A 2 /2 = 0.1 meV. In order to have strongly spin- 
polarizing QPC, the latter should support many trans- 
mitting channels at chemical potential [i ~ m*X 2 /2 S> 
Tilj. This condition can, equivalently, be written in terms 
of the width of the point contact Ay, see Eq. (f . as 
Ay ^> Ir. This is a realistic condition for typical ballis- 
tic constrictions. 

To conclude, we have proposed a way to polarize cur- 
rents in the ballistic regime by means of using electric 
gates to suppress transmission in the upper spin-orbit- 
split subband. The polarization is stronger when there 
are many transmitting channels in the lower subband. 
This is exactly the condition when the semiclassical ex- 
pansion in powers of h is applicable. An obvious advan- 
tage of our scheme is that we do not require the spatial 
modulation of the strength of spin-orbit interaction. Nei- 
ther do we need a restricted angle of incident electrons 
in order to have a polarized current. 
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